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Introduction 

The rare earth oxide-zirconia systems are used for a variety of applications, such as solid 
electrolytes and thermal barrier coatings. To gain a better insight into their behavior and 
further improve their performance, it is essential to understand the thermodynamic 
properties of these materials. These thermodynamic properties dictate the stability of the 
oxides in actual service conditions. Today thermodynamics involves both experimental 
measurements and development of self-consistent databases which contain the Gibbs 
energies of the individual phases (1). Parameters in various solution models are 
evaluated based on available experimental data. From a set of relatively simple 
equations, all phase boundaries and thermodynamic properties can be calculated for a 
particular system through minimization of the Gibbs energy under given conditions. 
Further, binary systems are building blocks for higher order systems, which may not be 
experimentally studied. 

Yttria-zirconia (Y2O3-Z1O2) is the most widely used of the rare earth oxide-zirconia 
systems. There are numerous experimental studies of the phase boundaries in this 
system. In this paper, we assess these data and derive parameters for the solution models 
in this system. There is current interest in other rare earth oxide-zirconia systems as well 
as systems with several rare earth oxides and zirconia, which may offer improved 
properties over the Y203-ZrC>2 system. For this reason, we also assess the ytterbia- 
zirconia (YF^Cb-ZrCh) and Y 2 03-Yb203-Zr02 system. Only limited experimental data 
are currently available for the Yb203-Zr02 binary and Y203-YI>203-Zr02 ternary systems. 

There are several calculations of the Y203-Zr02 system in the literature. Kaufman (2) 
approaches this type of psuedo-binary by treating each oxide unit (i.e. Y2O3 and Zr0 2 ) as 
a component. This readily allows development of quite good oxide databases with a 
minimum number of components. The experimental phase diagram for Y 2 0 3 -Zr02 is 
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dominated by a large cubic phase field for Zr 0 2 as well as a large cubic phase field for 
Y2O3. Degtyarev and Voronin ( 3 , 4 ) point out that the cubic Y2O3 phase can be derived 
from the ZrC>2 phase by removing two oxygen anions. Thus they treat the two cubic 
phases as one cubic phase with a miscibility gap. Du et al. ( 5 , 6 ) model the Y203-ZrC>2 
system using the miscibility gap simplification and treating YO1.5 and ZrC>2 as 
components. They do a thorough literature survey and derive their model coefficients 
using experimental data ( 7 ). They do not consider experimental data below —1300 K, due 
to likely difficulties in attaining equilibria. One of their assessments uses the liquid as the 
reference state ( 5 ) and the other uses the standard element reference state (6). Yokokawa 
et al. (8,9) calculate the Y 2 03-Zr0 2 diagram and phase diagrams for a number of other 
rare earth oxide-zirconia systems. They also treat the oxide units as components and 
derive solution parameters from both phase boundary data and correlations between ionic 
radii of metal oxides. 

No calculations were found for any ternaries or higher order systems. In this study we 
assess the Y2O3-Z1O2, Yb 2 03-Zr0 2 , and Y 2 03-Yb 2 03 binaries from available data. We 
re-assess the Y2O3-Z1O2 for two reasons — to build a database without the miscibility gap 
simplification and to include more low temperature data. We do not use the miscibility 
gap simplification for two reasons. First, the space groups of cubic Y2O3 and cubic Z1O2 
are not the same. Second, the use of separate phases makes convergence easier. The 
oxide units Y2O3, YF^Cb, and ZrC>2 are used as components. The model parameters are 
derived from experimental data. Then these three binary assessments are combined to 
form the ternary Y203-Yb203-Zr02 database. 

Review of Experimental Data 

Pure Components and Lattice Stabilities 

All data is in the form of Gibbs energies with the standard element reference state (SER). 
Transition heats and temperatures are given in Table 1 (a). These data are from the STGE 
database ( 10 ), which in turn were taken from IVTAN database ( 11 ). Adjustments are 
made to the transition temperatures of Y2O3 in order to maintain consistency with the 
Y2C>3-rich data of Stubican et al. ( 12 ), which is the most complete dataset for the Y2O3- 
rich side of the diagram. Further the IVTAN heat of transition between the two solid 
phases is 54 kJ/mol, which is rather high. A more reasonable estimate for a solid 
oxide/solid oxide transformation is 20 kJ/mol ( 13 ) and this is used in the database. 

The SGTE and IVTAN data for Yb 2 0 3 are incomplete. The SGTE database lists only 
one solid phase for Yb 2 03, whereas the IVTAN tables list two solid phases, but with no 
heat of transition between the two solid phases. In analogy with Y2O3, the lower 
temperature solid phase is taken as cubic and the high temperature solid phase as 
hexagonal and the heat of transition is estimated to be 20 kJ/mol. The transition 
temperatures are taken from IVTAN ( 11 ). 

Solution database development is based on the determination of lattice stabilities ( 1 ). For 
the pure elements, these have been determined and tabulated. However, when using 
oxide units as components, these must be estimated. Lattice stabilities are chosen to 
make the particular phase unstable over the entire range by adding -20 kJ/mol to the 
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stable phase’s Gibbs energy for a given temperature range. All lattice stabilities are listed 
in Table 11(b). 

Y 2 03 -ZrC >2 System — Phase Boundary Data 

There is a good deal of experimental phase boundary information on this system (12,14- 
38). The data through 1990 are reviewed by Du et al. (5,6). Figure 1(a) and (b) show 
two accepted experimental phase diagrams (14,15). Both diagrams are dominated by a 
large cubic Zr0 2 phase field. Both diagrams contain several invariant points, which are 
listed in Table II. The differences between the two experimental diagrams indicate the 
many controversies in this diagram. Du et al. (5,6) state that the difficulty in attaining 
low temperature equilibria make the low temperature eutectoids unreliable. Therefore 
they do not use them in their assessment. The Y 2 03 -Zr 0 2 system can be viewed as a 
rigid cation lattice with mobile anions. Due to the low diffusivity of the cations, phase 
changes are very slow, particularly at the lower temperatures. This fact must be 
considered in assessing experimental data. In the following discussion we examine 
available experimental data paying particular attention to the technique and reliability of 
the data. 

There are three measurements of the liquidus line (17-19). Skaggs (19) used a laser 
heating method to create a ‘pendant drop’. Temperature measurements were taken with a 
small amount of solid in the drop and then again with approximately half the drop solid. 
The average of these gave the liquidus temperature. However, these data showed a large 
amount of scatter and are not used in this assessment. Noguchi et al. (17) used a solar 
furnace in conjunction with a pyrometer. Cooling curves clearly revealed the solification 
point. Rouanet (18) has used a similar technique and his data are in basic agreement with 
Noguchi et al. (17). The data of Noguchi et al. (17) are thus used in this assessment. 

The Zr0 2 -rich side of the diagram is the area of the greatest practical interest. As noted, 
the equilibrium phases of Zr0 2 are monoclinic, tetragonal, and cubic as temperature 
increases. In addition, quenched samples also form a persistent, non-equilibrium 
tetragonal phase (32). The non-equilibrium tetragonal phase has a smaller c axis and 
generally higher Y 2 C >3 contents than the corresponding equilibrium tetragonal phase (37). 
Some of the discrepancies in reported experimental data are likely due to the formation of 
this metastable phase (32, 37). 

The most well-known transformation of Zr0 2 is the monoclinic to tetragonal 
transformation, which is accompanied by a 3-5% volume increase. Four measurements 
were found of this transformation (16, 24, 27, 31) — all show that the transformation 
temperature decreases with increasing Y 2 C >3 content. Duwez et al. (16) observed a sharp 
change in expansion on the tetragonal to monoclinic transformation. Pascual and Duran 
(24) used dilatometry and differential thermal analysis (DTA). Ruh et al. (27) used XRD 
and Yoshikawa and Suto (31) used dilatometry. In the mole fraction of Y 2 C >3 from 0.00 
to 0.02, these data all show reasonably good agreement. However, Duwez et al. (16) and 
Pascual and Duran (24) show the transformation temperature continues to decrease to a 
mole fraction of Y 2 C >3 - 0.05. Only the mole fraction of Y 2 03 from 0.00 to 0.02 data are 
used in this assessment. 
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There are the most data for cubic/tetragonal equilibria at the Zr 02 -rich side (12, 16, 21, 
23-25, 30, 36, 38). Duwez et al. (16) used XRD to determine the cubic/(cubic + 
tetragonal) boundaries at 1648 and 2273 K. Ruh et al. (27), Pascual and Duran (24), 
Scott (21) and Stubican et al. (12) have taken similar measurements with XRD. Ruhle et 
al. (25) and Lanteri et al. (23) used analytical electron microscopy to determine the 
tetragonal/(tetragonal + cubic)/cubic phase boundaries. Their data indicates these 
boundaries are temperature independent. Yoshikawa et al. (30) used electron microprobe 
analysis to determine these phase boundaries. The data points for all these investigations 
are collected in Figure 2 and show remarkably good agreement with each other at 
temperatures greater than -1600 K. At lower temperatures the mole fraction of Y 2 O 3 
differs by - 0 . 02 . 

The cubic phase field on the ZrCb-rich side occupies a large section of the diagram and 
extends to low temperatures (33). Stubican et al. (12) has proposed a low temperature 
eutectoid, as noted in Table II, which is included in this assessment. 

There is one confirmed intermediate compound in this system — Y 3 Zr 40 i 2 (12). This is 
an ordered phase, which decomposes to the disordered cubic solution. As shown in Table 
II, there is some controversy of its decomposition temperature. Other investigators have 
reported a pyrochlore-type compound of the form Zr 2 Y 207 (24). However the existence 
of this compound is controversial — many investigators have not observed it, despite 
careful studies to search for it (34). Therefore it is not considered in this assessment. 

As shown in the experimental phase diagram (Figure 1(a) and (b)), there are large cubic 
Zr 02 and cubic Y 2 O 3 regions, as well as a two phase region. Stubican et al. (12) have 
determined the phase boundaries of both sides of this two phase field by XRD. Pascual 
and Duran (24) have also used XRD to determine the cubic ZrC> 2 /(cubic ZrC >2 + cubic 
Y 2 O 3 ) boundary. Their data is a mole fraction of Y 2 O 3 -0.08 higher than that of Stubican 
et al. (12). Jayaratna et al. (36) have also measured these boundaries. Their data for the 
cubic Zr 02 /(cubic ZrC >2 + cubic Y2O3) boundary shows excellent agreement with that of 
Stubican et al. (12). The data of Stubican et al. (12) is used in this assessment. 

Finally, the Y 203 -rich side of the diagram is less well studied. It is dominated by the 
cubic Y 2 O 3 phase, which transforms to the hexagonal Y 2 O 3 phase at high temperatures. 
The liquidus data in this region and the phase boundary data is given by Stubican et al. 
(12). Again, this data is the most complete for this side. As noted, some corrections 
were made to the Y 2 O 3 lattice stabilities for consistency with these data. 

Y2C>3-Zr02 System — Thermodynamic Data 

There is only limited thermodynamic data on this system. Belov and Semenov (39) have 
measured the activity of Zr 02 and Y 2 O 3 at 2773 K using a Knudsen cell vapor pressure 
measurement technique. There are some inconsistencies between these data and the 
reported cubic/liquid phase boundary as shown in Figure 3. Assuming a pseudo-binary, 
there should be a constant activity in the (liquid + cubic) two phase region. There are two 
possible reasons for this inconsistency- — either the temperature measurement is off or 
since the system is not truly a binary, but rather a Y-Zr-0 ternary, P(C> 2 ) is varying in the 
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experiment. Following Du et al. (4), we only include the first three Zr0 2 -rich 
measurements in this assessment. 

One calorimetric measurement was found for the heat of formation of an intermediate 
compound of the formula Zr 2 Y 2 0 7 (40) and is given in Table III. However a compound 
of this stoichiometry has not been observed experimentally. First principles calculations 
are available for the heat of formation of ZrjY40i2 (41) and are also listed in Table III. 

Yb 2 03-ZrC>2 System — Phase Boundary Data 

There is much less experimental data available on this system and only phase boundary 
information was found. Only one liquidus measurement was found (18). The Zr0 2 -rich 
side of the at temperatures below the solidus line has been studied by Corman (26, 29). It 
is similar to the Zr0 2 -rich side of the Y 2 C>3-Zr0 2 phase diagram and also contains an 
intermediate compound Zr3Yb40i 2 , which disorders at a higher temperature than 
Zr3Y40i 2 as listed in Table II. 


Y203-Yb 2 03-Zr0 2 System 


The Y 2 03-Yb 2 03-Zr0 2 system appears to be the only pseudo-ternary system with Zr0 2 
and two rare-earth oxides, which has been experimentally studied. Corman (26) and 
Corman and Stubican (29) have examined this system at three temperatures from pure 
Zr0 2 to solutions with a mole fraction of Y 2 03 equal to 0.5 and a mole fraction of Yb 2 03 
equal to 0.5. Their diagrams at 1473, 1673, and 1923 K are shown in Figures 4(a), (b), 
(c), respectively. 

The ternaries are characterized by continuous tetragonal Zr0 2 , cubic Zr0 2 , and cubic 
(Y,Yb) 2 03 solutions. Further the intermediate line compound, Zr3(Y,Yb) 4 0i 2 , extends 
across the diagram when both ends are stable (Figure 4(a)). As temperature is raised 
above 1655 K, the Zr3Y40i 2 side disappears (Figure 4(b)). Above 1910 K, the compound 
disappears entirely (Figure 4(c)). 

Modeling 

As discussed, the oxide units are treated as components. The lattice stabilities are listed 
in Table 1(b). 

Each of the solutions are modeled with random mixing term and a Redlich-Kister 
polynomial (42) to describe the excess free energy: 

G„(*.T) = CD + RT 2>i >n x, + “G m 

i=I i=l 

=x,x 2 J j L(x,-x 2 ) j 

j=0 

Here x, is the mole fraction of component i, °G,(T) is the lattice stability of that 
component (Table 1(b)), 'L are the interaction parameters, R is the gas constant, and T is 


( 1 ) 

( 2 ) 
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the absolute temperature. In this assessment, only the first two terms of the Redlich- 
Kister polynomials were used: 

These interaction parameters were evaluated from the experimental data via the 
optimizing code Parrot in ThermoCalc* and also the optimizing code in WinPhad**. 
First the Y203-Zr02 diagram was obtained by fitting the liquidus and the approximate 
boundaries of the large cubic phase fields. The smaller phase regions were added by 
fine-tuning the parameters. A number of optimizations from different starting points 
were done to get the best fit with no artifacts. The Yb203-Zr02 diagram was done by 
using the parameters from the Y203-Zr02 system and making adjustments for the 
differences in the two systems. 

ThermoCalc uses the Muggianu approach (43) to extend the binaries to ternaries. This 
approach gives the correct limits when the ternary solutions reduce to the binaries. In 
addition a ternary interaction parameter was added of the form: 

ex ' en G m = x,x 2 x 3 0 L(tem) (4) 

These parameters were determined by ‘trial and error’ to get the best agreement with the 
experimental ternary sections. 

Results and Discussion 

The solution parameters for the Y203-Zr02 system and are shown in Table IV. In most 
cases the °L and *L parameters are needed; however in some cases only the °L parameter 
is needed. The complete calculated phase diagram is shown in Figure 5(a) with 
experimental data shown as symbols. Nearly all features of the experimental diagram 
(Figurel(a) and (b)) are reproduced including the five invariant points (Table II), the 
liquidus lines, and the cubic ZrC>2/(cubic ZrC>2 + cubic Y2C>3)/cubic Y2O3 phase 
boundaries. Only the tetragonal Zr02/(tetragonal Zr02 + cubic Zr02)/cubic Zr02 does 
not follow the experimentally measured phase boundaries, as shown in Figure 5(b). The 
highest temperature points are reproduced, but the calculated diagram does not reproduce 
the measured lack of temperature dependence of these boundaries. The appearance of the 
monoclinic solid solution at 1445 K also pushes these boundaries into higher Y2O3 
contents. 

Figure 6 shows the calculated activities in this pseudo-binary at 2773 K. Note the two 
phase (liquid + cubic Zr0 2 ) region with constant activity. The experimental data are 
shown as data points. The measured activity of Z1O2 was used in the assessment and 
shows very good agreement. The activity of Y2O3 was not used and shows approximate 
agreement with the calculations. 

The enthalpy of formation for Zr3Y40i2 is also extracted from this assessment from: 

3 Zr0 2 + 2 Y2O3 = Zr3Y 4 0 I2 (4) 


’ ThermoCalc Ver. N, ThermoCalc Software, Stockholm, Sweden 
*' WinPhad Pro 2.0, Computherm LLC, Madison, WI. 
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The value from the assessment (- 45.9 kJ/mol) shows order of magnitude agreement with 
the calculations from first principles (Table III, ref. ( 41 )). 

The parameters for the Yb203-Zr02 system are given in Table V and the calculated phase 
diagram is given in Figure 7 . There is good agreement with the limited experimental 
data. Due to these limited data, the monoclinic ZrC>2 phase was included only as a line 
compound. The calculated tetragonal ZrC>2/( tetragonal ZrC>2 + cubic Zr02)/cubic ZrC>2 
shows better agreement with experimental data than in the Y2C>3-Zr02 system, as shown 
in Figures 7 and 5 , respectively. 

The third binary needed for this study is the Y203-Yb203 system. This was treated as an 
ideal solution, in accordance with the limited experimental data ( 44 ), which shows a 
continuous cubic solid solution at 1923 K. 

Next the binaries were combined to form the Y203-Yb2C>3-Zr02 psuedo-temary. 
Additional ternary interaction coefficients were included and given in Table VI. 
Calculated sections for 1473 , 1673 , and 1923 K are shown in Figures 8 (a), (b), and (c), 
respectively. The intermediate line compound Zr3(Y,Yb)40i2 shows the experimentally 
observed behavior (Figures 3 (a-c)). The tetragonal ZrC>2, cubic Zr 02 , and cubic 
(Y,Yb)2C>3 phases also show continuous solid solutions, as observed experimentally. 

Conclusions 

Experimental data for the Y 2 C> 3 -Zr 02 , Yb203-ZrC>2, and Y203-Yb203-Zr02 systems has 
been examined. For the Y 2 0 3 -Zr 0 2 psuedo-binary a good deal of phase boundary data 
are available, but only limited thermodynamic data. For the other systems, only limited 
phase boundary data is available. A Calphad-type database was developed for both 
pseudo-binaries and put together for the pseudo-temary. 
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Phase Transition 

Temperature of transition (K) 

An-arnH (kJ/lT10l) 

Reference 

Zr0 2 monoclinic — > tetragonal 

1445 

8400 

(ip 

Zr0 2 tetragonal — > cubic 

2620 

21400 

(id 

Zr0 2 cubic liquid 

2983 

111400 

(ip 

Y 2 O 3 cubic — ► hexagonal 

2550 

54000 

(ip 


2599 

20000 

J12! 

Y 2 0 3 hexagonal ->■ liquid 

2712 

115000 

(ii) 


2782 

114000 

J1D1 

Yb 2 03 solid — ► liquid 

2789 

130000 

(ii) 

Yb 2 0 3 cubic -> hexagonal 

2663 

20000 

(ii)** 

Yb 2 0 3 hexagonal -> liquid 

2703 

12000 

( 11 )** 


•Values from reference 1 1 modified to give agreement with data of reference 12. 
••Transition temperatures from reference 1 1 and heats are estimated. 


Table 1(a). Temperatures and enthalpies of transitions to define lattice stabilities. 


Phase 

Abbrev 

iation 

Temperature 
Range (K) 

Gibbs energy (J/mol) 

Zr0 2 — monoclinic 

mz 

298-1445 

-1 125582.53+417.508432*T-68.32783*T*LN(T)- 
.00454 1 5885 *T**2+l. 72444 167E-10*T**3+671 863. 5 *T**(-1) 



1445-2620 

1129291 49+484.935707T-78. 1 *T*LN(T)-1 .0845475E-1 6*T**2 
+7.8 1505 167E-21*T**3-7.622945E-08*T**(-1) 



2620-2983 

- 1 1 34269.49+50 1 ,790474*T-80*T*LN(T)-2.0 126725E- 
14*T**2+9.19196667E-19*T**3-5.127845E-05*T**(-1) 



2983-6000 

- 1 1 93929.49+68 1 .8041 69*T- 1 00*T*LN(T)+3 .041038E- 
16*T**2-8.75008833E-21*T**3+2.46697E-06*T**(-1) 

Zr0 2 — tetragonal 

tz 

298-6000 

GnttfT) + 8400-5.8 13 14879*T 

Zr0 2 — cubic 

cz 

298-6000 

G ml (T)+2 1 400- 1 0.7749809*T 

Zr0 2 — liquid 

Lz 

298-6000 

G ro (T)+lt 1400-40.9459497*T 

Zr02 — cubic Y 2 0 3 phase 

czy 

298-6000 

G mz (T)+21500 

Z 1 O 2 — hexagonal 

hz 

298-6000 

Gmz(T)+22500 

Y 2 0 3 — cubic 

cy 

298-1100 

2.98140E+02 -1945493.95+680.542373*T-1 13.9721*T*LN(T) 

-.0 1276 1235*T**2+1. 49941 783E-06*T**3+813690.5*T**(-1) 



1100-2100 

-2012427.83+1275.7521 1*T-98.0735*T*LN(T)+.03498891*T**2 
-3.688465E-06*T**3+10928115*T**(-1) 



2100-2599 

-2274854.5 1+2766. 12583*T-94.0992*T*LN(T)+.l 003 1335*T**2 
-7.7830 1833E-06*T**3+77093000*T**(-1) 



2599-2782 

-201 16 17.72+ 1043.70289*T-160*T*LN(T)-5.08929E-13*T**2 
+2.4271 9833E-17*T**3-.001 1431 1*T**(-1) 



2782-6000 

-2 1 20097.72+ 1 399.92055*T-200*T*LN(T)+4.6124375E- 
16*T**2-1.38251267E-20*T**3+3.154309E-06*T**(-1) 

Y 2 0 3 — hexagonal 

hy 

298-6000 

G cv (T)+20000-7.69527*T 

Y 2 0 3 — liquid 

Ly 

298-6000 

G CV (T)+ 1 34000-48. 67307*T 

Y 2 Oj — monoclinic 

my 

298-6000 

G CV (T)+22100 

Y 2 0 3 — tetragonal 

ty 

298-6000 

G cv (T)+2t 100 

Y 2 0 3 — cubic Zr0 2 phase 

! cyz 

298-6000 

G CV (T)+20100 

Yb 2 0 3 — cubic 


| 298-2663 

! 

- 1 858343 .67+747 .94 1586*T- 1 30. 1 85 1 *T*LN(T) 

-.001 8969 155"'T**2-7.66497333E-12*T**3+709650*T**(-1) 



2663-2798 

- 1 923755 .38+1 002.70962*T- 1 60*T*LN(T) 



i 2798-6000 

-2035675 .3 8+ 1 360. 1 7602*T-200*T*LN(T) 

Yb 2 0 3 — hexagonal 

n 

298-6000 

G CV b(T)+ 20000-7.5 10327*T 

Yb 2 0 3 — liquid 

n 

298-6000 

G cvb (T)+ 140000-5 1 .823473 *T 

Yb 2 0 3 — monoclinic 


298-6000 

G cvb (T)+ 22100 

Yb 2 0 3 — tetragonal 


298-6000 

G CV b(TH 21100 

Yb 2 0 3 — cubic Zr0 2 phase 



298-6000 

G CV i,(T)+20 1 00 


Table 1(b). Lattice stabilities. 
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Invariant 

Temperature (K) 

Composition x(Y 2 03 ) 

Reference 

Comments 

Eutectic 
liq = cz + cy 

2688 

0.817 

( 12 ) 



2603 

0.76 

jm I 



2646 

0.825 

(19) 


Peritectic 
cy = liq + hy 

2720 

0.778 

( 12 ) 

Introduced, not measured 

Euctectoid 
cz = delta y + cy 

1650 

0.435 

(27) 


Euctectoid 
tz = mz + cz 

838 

0.039 

( 21 ) 

Difficult to attain equilibria at 
this low temperature 

Euctectoi 

cz = mz + delta y 

<673 

0.20-0.30 

( 12 ) 

Difficult to attain equilibria at 
this low temperature 
Approximate euctectoid 
discussed in ( 12 ) 

Decomposition 
Temperature of delta y 

1523 

0.4 

( 12 ) 



1655 

0.4 

(35) 



1648 

0.4 

(24) 


Eutectic 
liq = cz + cyb 

2691 

0.85 

(19) 


Decomposition of 
delta yb 

1898 

0.4 

(29) 

- 



Table II. Experimentally determined invariant points in the Y 203 -Zr (>2 and Yb 2 03 -Zr 02 
systems and decomposition temperatures of delta phase. 


Compound 

Temperature 

(K) 

Method 

A(H°(from oxides) — 
kj/mol 

ArH°(from elements) — 
kJ/mol 

Reference 

Y 2 Zr 2 07 

298 

Calorimetry 

-47.7 

-4185.4 

(36) 

Zr3Y40i2 

0 

First Principles — 
Generalized Gradient 
Approximation 

-20.9 

-6454 

(37) 

Z 13 Y 4 O 12 

0 

First Principles — 

Local-Density 

Approximation 

-14.7 

-6448 

(37) 

Table III. Entl 

halpies of 1 

'ormation of intermediate compound. 


Phase 

Parameter 

Value (J/mol) 

Liquid 

°L 

98445 -45.67T 


'L 

95521 -30.52T 

Cubic Zr0 2 

°L 

-144364 + 24.60T 


'L 

-1 80000 + 62.25T 

Tetragonal Zr0 2 

°L 

-22324 - 32.95T 


'L 

-240000 + 50T 

Monoclinic Zr0 2 

°L 

-90000 - T 

Hexagonal Y 2 C >3 


-18000 

Cubic Y 2 0 3 

°L 

-30735 -20.97T 


'L 

99126 -66.03T 

Zr3Y40i 2 

A+BT 

-45966 + I1.28T 


Table IV. Interaction parameters for Y 2 C>3-Zr02 system. 
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Phase 

Parameter 

Value (J/mol) 

Liquid 

°L 

18466- 15.72T 


l L 

-42812 + 26.53T 

Cubic Zr0 2 

“l 

-116088+ 12.67T 


‘L 

-296755 + 98.69T 

Tetragonal ZrC >2 

°L 

-22326+ 33. 00T 


'L 

-240024 + 50T 

Monoclinic Zr0 2 

°L 

-90000 - T 

Hexagonal Yb 2 03 

°L 

-18000 

Cubic Yb 2 03 

°L 

-30735 - 20.97T 


'L 

125000 -66.03T 

Zr 3 Yb40 , 2 

A+ BT 

-50064 + 1 IT 


Table V. Interaction parameters for Yb^C^-ZrCh system. 


Ternary Phase 

L(Y 2 0 3 -Yb 2 0 3 -Zr0 2 ) (J/mol) 

Tetragonal Zr0 2 

0 

Cubic Zr0 2 

110000 

Zr 3 (Y,Yb) 4 Oi 2 

-8000 

Cubic (Y,Yb),0 3 

-100000 


Table VI. Ternary interaction parameters for Y 2 C> 3 -Yb 203 -Zr 02 
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Figure 2 — Experimental data for the tetragonal/(tetragonal + cubic)/cubic 
phase boundaries. 
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Figure 3. Activity of Zr0 2 and Y 2 0 3 at 2773 K from Belov and Semenov (35) 
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# Noguchi 
Q Stubican (1978) 
x Ruhle 

# Srivastava 

A Stubican (1984) 

# Yoshikawa 
t Pascual 



MOLE FRACTION 03Y2 


9 Noguchi 
o Stubican (1978) 
x Ruhle 

# Srivastava 

* Yoshikawa 


Figure 5. Calculated phase diagram for the Y 2 0 3 -Zr0 2 system, (a) Entire 
diagram (b) Enlargement of Zr0 2 -rich region. 
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Figure 6. Plot of calculated activity in Y 2 0 3 -Zr0 2 system vs composition, 
showing experimental measurements (35). 



# Rouanet 
□ Corman 


Figure 7. Calculated Yb 2 0 3 -Zr0 2 phase diagram. 
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Figure 8. Calculated isothermal sections for Y 2 C> 3 -Yb 203 -Zr 02 at (a) 1473 (b) 1673 and 
(c) 1923 K. 



